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ON THE ABUNDANCE OF COLLAPSED OBJECTS 
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ABSTRACT 

The redshift dependence of the abundance of coUapsed objects places strong constraints on cosmo- 
logical models of structure formation. We apply a recently proposed model describing the anisotropic 
collapse of inhomogeneous spatial domains. Compared with the spherical top-hat model, this generic 
model leads to significantly more collapsed objects at high redshifts: at redshift one and on the scale of 
rich clusters a factor of 65. Furthermore, for a fixed normalization of the initial fluctuation spectrum 
(cTg = 1), we predict four times as much presently collapsed objects on the mass-scale of rich clusters 
within the standard CDM cosmogony, compared to the spherical collapse. 

Subject headings: galaxies: clusters: general — galaxies: abundances 
large-scale structure of universe 



■ cosmology: theory 



1. INTRODUCTION 

The abundance of galaxy clusters has been studied by 
Press & Schechter (1974) based on the spherical self- 
similar collapse. If this spherical domain is initially over- 
dense, one is led to a model of collapse, featuring a scale- 
factor (proportional to the radius of the domain) that, 
in the course of evolution first increases with the Hubble 
flow, then attains a vanishing time-derivative interpreted 
as the "decoupling from the global Hubble-flow", and fi- 
nally shrinks to zero interpreted as "the collapse of the 
spherical region" . On a given spatial scale the spherical 
top-hat model consistently describes an embedded "small 
Friedmann model" , its radius (scale-factor) also obeying 
Friedmann's equations but for a mass different from the 
mass given by the background density. This top-hat model 
is a sub-case of the general expansion law discussed below. 
Deviations from the spherical collapse in a realistic situa- 
tion have been studied in numerous works (some of which 
shall be cited later on). Notwithstanding, results from 
these works have not led to a simple alternative model for 
generic anisotropic collapse of inhomogeneous spatial do- 
mains, and are mostly considered as (local) refinements 
of "rough but sufficiently accurate" top-hat estimations. 
In this situation we have recently proposed a model for a 
generic collapse as a spin-off from studies of the "backreac- 
tion problem" in Newtonian cosmology. We shall demon- 
strate that on large mass-scales this generic model leads 
to substantially different results that cannot be regarded 
as mere refinements of the top-hat model. 

The generic collapse model is based on a general ex- 
pansion law in Newtonian cosmology for the dust model, 
obtained by averaging Raychaudhuri's equation on mass- 
conserving spatial domains, as summarized below. Let us 
denote with the spatial averaging in Eulerian space, 
e.g., for a spatial tensor field yl(x, t) — {Aij (x, t)} we sim- 
ply have the Euclidean volume integral normalized by the 
volume of the domain: {A)-jy = 1/V{t) j^d^xA^X-.t). Es- 
pecially we are interested in the time evolution of the vol- 



ume V — |P| of the domain T), modeled by the domain- 

1/3 

dependent scale-factor aT>{t) — {V{t)/V{ti)) ' . By av- 
eraging Raychaudhuri's equation for mass-conserving do- 
mains, i.e. domains identified in Lagrangian space, Buchert 
& Ehlers (1997) obtained a general expansion law describ- 
ing the evolution of the volume of a domain, via the evo- 
lution of the scale-factor ax>: 



2,— +A^G{q)^-K^Qi 



(1) 



with Newton's gravitational constant G, the cosmological 
constant A, and the "backreaction term" Q-p, which fea- 
tures positive definite fiuctuation terms in the parts of the 
velocity gradient: 

Qv^l{{e-{e)pr)^ + 2{c.'-a')^, (2) 

with the expansion rate 9, the rate of vorticity lo, and the 
rate of shear a. For = this equation equals one of the 
Friedmann equations for the scale-factor a{t) — ax>{t) in 
a homogeneous and isotropic universe with uniform den- 
sity qh = {q)j}- The key-difference of the evolution of 
generic domains compared with that described by the top- 
hat model may be summarized as follows: the collapse of 
a generic domain is triggered not only by an over-density, 
{^)v ~ {{b)v ~ Qh)Iqh, but also by fluctuations in the 
velocity gradient encoded in the "backreaction term" Qp, 
most prominently by the averaged expansion and shear 
fluctuations. Consistently, Q-p vanishes for spherically 
symmetric flows inside spheres. 

The part still needed to solve the general expansion 
law (1) is a model for the time evolution of Qp. Buchert 
et al. (2000) (BKS) calculated Qv{t) for the spherical and 
plane collapse, and also provided the results based on the 
growing mode solution of the Eulerian linear approxima- 
tion and on the Lagrangian linear approximation restricted 
to Zel'dovich's approximation. The approximate Qvit) 
based on the "Zel'dovich approximation" is appropriate, if 



one wants to trace generic initial conditions into the mildly 
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nonlinear regime. It is a powerful property of the general 
expansion law (1) combined with this backreaction model 
that the plane-symmetric collapse and the spherical col- 
lapse are included as exact sub-cases. Important for our 
work is that the "backreaction term" ^^(I/;, IIj^, 111^, t) 
now only depends on known time-dependent functions, 
and the invariants Ij^, Il/f , TLIr of the initial velocity gra- 
dient Vij averaged in an initially spherical volume with 
radius R (see BKS). For a Gaussian random field the 
volume-averaged invariants In, II/j, Illij are uncorrelated 
and the statistical (ensemble) average of each of them is 
equal to zero for any domain (BKS). But the fluctuations 
of, e.g., cr|(i?) — E are non-zero and may be cal- 
culated from the power spectrum. Specifically, for the 
volume-averaged first invariant this is straightforward: 

aUR) - [ d^k Pik)WR{k)^. (3) 

P{k) is the initial power spectrum and W}i{k) is the 
Fourier transform of the top-hat window with radius R. 
Similarly, an (R) and am (R) may be related to the power 
spectrum of the density fluctuations (see BKS). The vari- 
ance crj(i?) is equal to the well-known mean square fluc- 
tuations of the initial density contrast field. 

Table 1 

The r.m.s. fluctuations o-i(-R), o"ii(-R), and criii(R) for an 

INITIAL DOMAIN OF RADIUS -R ARE GIVEN: THE CALCULATION WAS 
BASED ON A STANDARD CDM POWER SPECTRUM. TO MAKE THESE 
NUMBERS MORE ACCESSIBLE, ALSO THE LINEARLY EXTRAPOLATED 
MEAN FLUCTUATIONS a{to)c7i(R) FOR A DOMAIN WITH SCALED 
RADIUS a{to)R AT PRESENT TIME ARE GIVEN. M IS THE TOTAL 
MASS INSIDE SUCH A DOMAIN. 



a{tQ)R [/i-^Mpc] 


5 


8 


12.5 


a{to)ai 


1.5 


1.0 


0.6 


R [/i-iMpc] 


0.025 


0.040 


0.062 


0.15 


0.6 


2.3 


CTI XlO^ 


7.5 


5.0 


3.0 


an xlO^ 


35 


16 


7.4 


am xlO^ 


15 


4.2 


1.0 



In the following we assume that the background evo- 
lution follows an Einstein-de Sitter model with h = 
0.5 = ffo/lOO/i^ikm/s/Mpc. The scale-factor of the 
background model is normalized to unity at the initial 
time, hence a{ti) = 1 and a{to) — Zi + 1, and the ini- 
tial density fluctuations are modeled as a Gaussian ran- 
dom field with a Cold Dark Matter (CDM) power spec- 
trum (Bardeen et al. 1986). We employ the normalization 
a{to)ai{8h~^Mpc/a{to)) — as = I- For such a situation 
the fluctuations cti, cth, ctih are given in Table 1 for sev- 
eral radii of the initial domain. The normalization of the 
power spectrum enters in a linear, quadratic, and cubic 
way into ai{R), au{R) and am{R), respectively. 

2. GENERIC COLLAPSE 

As an operational definition for the redshift of collapse 
Zc we use the redshift when ax>{zc) approaches zero with 
ciT>{zd) < 0. Our interest is whether a domain did al- 
ready collapse until present. Clearly, this will depend on 
the initial conditions and on the radius R of the initial 
spherical domain, quantified by the r.m.s. fluctuations of 



the volume-averaged initial invariants f7i(i?), aii{R) and 
o'iii(i?). To calculate the distribution of Zc we assume 
that for a fixed radius R the volume-averaged initial in- 
variants lij, Ur and TLLr are Gaussian random variables 
with zero mean and variance ai{R), aii{R) and am{R), 
respectively. This is correct for Ir, but certainly an ap- 
proximation for Ilfl and TLLr (see the discussion below). 
We choose values for Ir, IIr and 111^, according to these 
probability laws. Then we solve the general expansion 
law (1) numerically for a-D as outlined in BKS starting 
at a redshift of Zi = 200 and calculate the redshift Zc, if 
the domain did collapsed up until z = 0. We iterate this 
procedure to determine the distribution of Zc- Starting 
at Zi = 1000 or 50 with accordingly rescaled fluctuations 
leads to nearly identical results. 

For either spherically symmetric {Qv — 0), plane- 
symmetric (see Eq. (44) in BKS), or generic initial con- 
ditions (see Eq. (42) in BKS), we estimate the probability 
p{zc)dzc that a domain collapses at Zc, using the Monte- 
Carlo procedure outlined above. Hence, dzc p{zc) is 
the fraction of domains which collapsed until present, and 
F(zc) = Jj'^ dz p{z) is the fraction of the domains which 
already collapsed up to a redshift of Zc- 
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Fig. 1. — The fraction F{zc)/F{0) of collapsed regions with co- 
moving "radius" of 8/i~^Mpc is shown: generic initial conditions 
(solid line), the plane collapse (dashed line), the spherical collapse 
(dotted line), and the generic collapse for restricted initial condi- 
tions with Ijj < —1.69 (dashed dotted line). In the inset plot p(zc) 
is shown. 

In Fig. 1 both p{zc) and F{zc) / F{0) are shown for a do- 
main with initial radius of 0.04/i~^Mpc, a comoving "ra- 
dius" of 8ft.~^Mpc corresponding to a mass-scale of rich 
clusters. The most striking feature is that the generic 
model leads to significantly more collapsed regions at high 
redshift in comparison to the spherical model. Moreover, 
only 5% of the domains collapsed in the spherical model, 
whereas in the plane-symmetric case 16%, and in the 
generic case already 20% of the domains collapsed. Rich- 
stone et al. (1992) studied F{zc)/ F{0) using the spherical 
model (see Bartelmann et al. 1993 for consistent initial 
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conditions). Considering domains with an initial radius 
5h~^Mpc/a{to) wc find the same result as Bartelmann 
et al. (f993) for the F{zc)/F{0) in the spherical model. 
Initially over-dense domains with —I/? = (5)^ > 1.69 will 
always collapse in the top hat model. In an Einstcin- 
de Sitter universe the top-hat model predicts that domains 
collapsed only recently (Fig. 1). As already mentioned, 
Ilfl and Illfl are stochastically independent of Ir for a 
Gaussian random field (see BKS) and we may investigate 
the generic case undcir the condition that ((5)^ > 1.G9. 
Now the majority of these initially over-dense domains 
already collapsed before a redshift of one. 

To investigate how strong our results depend on the as- 
sumption of a Gaussian law for aii{R) and (Jiii(i?) we 
repeat our calculations keeping ai{R) fixed but doubling 
(halving) CTii(i?) and crui{R)- The increased fluctuations 
mimick the effect of a pronounced tail of the distribution 
and consequently lead to even stronger deviations from 
the top-hat model. For reduced fluctuations the results 
stay in-between the generic model and the plane collapse. 
In both cases the redshift dependence of the normalized 
F{zc)/F{0) shows nearly no difference compared to the 
generic case given in Fig 1. 

Up to now we focused on the collapse of domains with an 
initial radius R = 8/i~^Mpc/ {zi + 1) containing the mass 
0.6 X 10^^ h~^MQ. In the following we will show how this 
accelerated collapse depends on the mass-scale. Using our 
Monte-Carlo procedure outlined above, wc calculate the 
fraction i^generic(^c)/-Fsphericai(-Zc) for a given mass-scale 
quantifying the deviation of the generic from the spherical 
collapse. 
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Fig. 2. — The ratio ^gen(2c)/^^Bph(2c) is plotted against the red- 
shift of collapse Zc for domains with a comoving size of 5^~^Mpc 
(solid line), 8fe~^Mpc (dashed line), and 12.5/i~^Mpc (dotted line). 

In Fig. 2 we see that at redshift z = the spheri- 
cal model under-predicts the abundance of galaxy clus- 
ters on the mass-scales considered. For smaller mass- 
scales (< 10^^h~^MQ) this discrepancy gradually van- 
ishes. However at higher redshifts and for larger mass- 
scales we predict significantly more collapsed domains, e.g. 



a factor of 65 at z = 1 and at a mass-scale of rich clus- 
ters (6 X lO^'*/i~^M0). We want to emphasize that the 
"Zel'dovich approximation" entering our calculations in 
the generic case becomes increasingly more accurate for 
larger domains. As discussed by Kitayama & Suto (1996) 
the SCDM model with as = I over-predicts the present 
abundance of galaxy clusters. Our main interest is in 
the comparison of the spherical collapse model with the 
generic anisotropic collapse model, using SCDM as a well 
studied reference model. 

3. DISCUSSION 

We examined the collapse of initially spherical domains 
embedded into a generic inhomogcneous cosmology. We 
combined the general expansion law for such domains with 
a "backreaction model" based on the "Zel'dovich approx- 
imation" to obtain explicit solutions for the scale-factor 
of such domains. Although we approximated the back- 
reaction term, the resulting model is more general than 
the volume deformation based on the averaged "Zel'dovich 
approximation". It still contains the spherical top-hat 
model as exact sub-case and is also exact in the orthog- 
onal case of plane-symmetric collapse. Assuming Gaus- 
sian initial conditions of a standard Cold Dark Matter 
cosmogony we calculated the abundance of collapsed do- 
mains. For a normalization of the initial fluctuation spec- 
trum with (78 = 1, collapsed domains on the mass scale 
of rich clusters are more abundant by a factor of four. 
Clearly, these results depend on the adopted normaliza- 
tion and background model. However, nearly independent 
of the normalization is the disagreement between the top- 
hat model's prediction and that of the generic anisotropic 
collapse model concerning the abundance of collapsed ob- 
jects at high redshift: we obtained 65 times more collapsed 
objects at a redshift of one on the scale of rich clusters. 
Compared to the spherical model, which focuses only on 
the mass over-density, our generic collapse model consid- 
ers also the fluctuations in the expansion rate as well as 
the effect of shear fields on the collapse (see Eq. (2)). Our 
results indicate that on large mass-scales the shear fields, 
caused by internal and external mass concentrations, ac- 
celerate the collapse. 

As a rule, our model mimicks that of N-body results 
at least down to the scale where the local (truncated) La- 
grangian approximation reproduces the results of N body 
runs (Weifi ct al. 1996). Only recently conducted large 
N— body simulations were able to provide enough dynamic 
range for the determination of the abundance of collapsed 
domains on the high mass-end (Jenkins et al. 2001; Bode 
et al. 2001). Similar to our results both studies show 
that the spherical model under -predicts the abundance 
of collapsed objects at high mass-scales. The abundance 
of high-mass objects drops significantly with increasing 
redshift. Hence the statistical significance of these results 
from N-body simulation is limited by the small number 
of massive halos. Nonetheless Bode et al. (2001) observe 
a similar growing deviation from the predictions of the 
spherical model with increasing redshift. Clearly, the pre- 
diction of an increase in the abundance of collapsed ob- 
jects by a factor of 65 at ^; = 1 needs to be verified against 
future simulations. If we normalize the model to the ob- 
served cluster abundance with as = 0.6 (Viana & Liddle 
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1996), the relative difference between the spherical and the 
generic collapse becomes even stronger (see also Governato 
et al. 1999). 

Sheth et al. (2001) estimated the mass function assum- 
ing an anisotropic ellipsoidal collapse, where the time evo- 
lution of the half-axes is determined by the "Zel'dovich 
approximation" . Our model may be viewed as a general- 
ization of this ansatz, where the domain is not restricted 
to be ellipsoidal throughout its evolution. Several other 
approaches are also based on the "Zel'dovich approxima- 
tion" , sometimes the third-order Lagrangian perturbation 
approximation is employed. In these "local" approaches 
the collapse of a domain is associated with a diverging lo- 
cal matter density (see e.g. Bartelmann et al. 1993, Monaco 
1997a,b, Lee & Shandarin 1998, and Engineer et al. 2000). 
The redshift of collapse is calculated from the distribution 
of the eigenvalues of the deformation tensor (Doroshkevich 
1970). The eigenvalues of the deformation tensor are de- 
termined locally at each point in space from a smoothed 
density field, where the smoothing scale determines the 
mass-scale. Compared to this pointwise treatment, our 
model is based on the time evolution of a finite domain. 
The invariants, averaged over the volume of the domain, 
determine the time-scales of the collapse. Therefore, our 
approach explicitly takes the spatial correlations of the 
initial velocity and density fields into account. 

We calculated the fraction of collapsed domains for a 
given mass-scale, regardless whether these domains are 
included in a larger collapsed region. This is the cloud- 
in-cloud problem as tackled by Bond et al. (1991) and 
Jedamzik (1995) for the spherical collapse (see also Yano 
et al. (1996) taking into account spatial correlations). Due 
to the drop of the mass-function at the large-mass end, the 
cloud-in-cloud problem mainly affects the smaller mass- 
scales. Lee & Shandarin (1998) showed in their calcula- 
tions of the mass-function that considering the cloud in 
cloud problem mainly changes the normalization, but not 



the shape. Hence, for the mass-scales considered in this 
Letter we do not expect a significant impact on the time- 
scales of the collapse. 

Still, the question remains, how the abundance of col- 
lapsed domains in this generic model may be related to 
the abundance of observed clusters. In the spherical model 
one can use the virial theorem to estimate the over-density 
where the collapse stops and the system "virializes" . How- 
ever, for small generic domains embedded into a fluctu- 
ating background and admitting non negligible generic 
shear fields on the size of the domain, neither the spheri- 
cal assumption nor the virial theorem for isolated systems 
(Chandrasckhar & Lee 1968) holds. 

Both, the spherical and the generic dust model lack 
forces that could compensate the collapse. Besides bary- 
onic physics (Gunn & Gott HI 1972) it would be manda- 
tory to include multistream forces opposing the collapse 
(Buchert & Dommguez 1998). While a spherical collapse 
is never realized, the generic model includes the effect of 
shear fields and fluctuations of the environment of the do- 
main to describe a realistic collapse situation. 

The results presented in this Letter have an obvious 
and important implication for ongoing high-redshift clus- 
ter survey projects using large telescopes such as Subaru 
and VLT. Our study shows the quantitative importance 
of modeling a realistic collapse situation for an accurate 
theoretical prediction of the mass-function. 
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